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Abstract 

Electron cooling that results when a bunch of electrons overlaps a 
bunch of ions , with both bunches moving at the same velocity, may 
be considered to be an intrabeam scattering process. The process is 
similar to the usual intrabeam scattering, where the ions scatter from 
each other and usually results in beam growth. An important differ- 
ence is that in electron cooling the mass of the ion is different from 
and much larger than the mass of the electron. This difference con- 
siderably complicates the intrabeam scattering theory. It introduces 
a new term in the emittance growth rate, which vanishes when the 
particles are identical and their masses are equal, and can give rise to 
emittance cooling of the heavier particles . The term that gives rise to 
beam growth for the usual intrabeam scattering is also present but is 
much smaller than the cooling term when one particle is much heavier 
than the other. This paper derives the results found for the emittance 
cooling rates due to the scattering of the ions in the ion bunch by the 
electons in the electron bunch. 

1 Introduction 

Electron cooling that results when a bunch of electrons overlaps a bunch of 
ions , with both bunches moving at the same velocity, may be considered to 
be an intrabeam scattering process. The process is similar to the usual intra- 
beam scattering, Ref.[l] where the ions scatter from each other and usually 
results in beam growth. An important difference is that in electron cooling 
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the mass of the ion is different from and much larger than the mass of the 
electron. This difference considerably complicates the intrabeam scattering 
theory. It introduces a new term in the emittance growth rate, which van- 
ishes when the particles are identical and their masses are equal, and can 
give rise to emittance coohng of the heavier particles . The term that gives 
rise to beam growth for the usual intrabeam scattering is also present but is 
much smaller than the cooling term when one particle is much heavier than 
the other. 

This paper derives the results found for the emittance cooling rates due 
to the scattering of the ions in the ion bunch by the electons in the electron 
bunch. The derivations given below makes considerable use of the results 
found in two previous papers, Ref.[2] and Ref.[3] 



2 The f{x,p) distribution and the scattering 
rate 6N 

The ions are contained within a bunch and their distibution is given by 
fa{xa,Pa) where Nafa{xa,Pa) IS the number of ions in d?Xad?Pa- Na is the 
number of ions in the bunch. 

j Sxad^Pa fa{Xa:Pa) = 1 

The distribution of the electrons in the electon bunch is given by fb{xb,pb) 
and Ni, is the number of electrons in the electron bunch. Let SNa be the 
number of ions with momentum, in (Ppa and space coordinate x in (Px 
which are scattered by the electrons with momentum pb in d^pb which are 
also in d^x, in the time interval dt , into the solid angle dfl' corresponding to 
the direction p'^. Then 5Na is given by, Ref.[2], 

NaNbaabdn'^^^^fa{x,Pa)fb{x,Pb)F{pa,Pb)d^xdt 

la lb 

[{PaPbf - mlmlY/^ 
rriamb 

Gab is the scattering cross section for the scattering of the ions from the 
electrons. In the expression for F{pa,Pb), we have put c — 1. F{pa,Pb) has 
the dimensions of a velocity. 



F{Pa,Pb) = 
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For completeness sake this result is given in the form which is valid in any 
CS. For the electron cooling problem for RHIC, one can do all the calclations 
in the Rest CS, which is the CS moving along with the two bunches. In the 
Rest CS, the central particle in either bunch is at rest and the motion of the 
motion of the particles may be treated non-reletavistically.In the Rest CS , 
one may put 7a = 76 = 1 and 

F{Pa,Pb) = \v*a - Vb\ 

3 Growth rates for < PiaPja > 

Following Bjorken and Mtingwa, Ref.[4], cooling rates will first be given for 

< PiaPja >• where the <> indicate an average over all the particles in 
the bunch. From these one can compute the growth rates for the average 
emittances of the ions, < e^a >. In a scattering event, where an ion with 
momentum Pa scatters off an electron with momentum p^, the momenta will 
change to p'^ and p'^. Let 5pia represent the change in pi^ in the collision, and 
similarly for S(piaPja)- Then 

Spia = P'ia-Pia 
SiPiaPja) = P'iaP'ja - PiaPja (2) 

Using the scattering rate given by Eq.(l), one can now compute 5 < 
PiaPja > in the Rest CS, 

S<{piaPja)> ^ ^bj d^xd^Pad^Pbfa{x,Pa)fb{x,Pb)\Va-Vb\ 

(TabdO! 5(piaPja) dt 
SiPiaPja) = {p'iaP'ja - PiaPja) (3) 

The 11- dimensional integral in Eq.3 can be reduced to a 3-dimenional 
integral for gaussian distributions, if one notes that in the Rest CS 
depends on iTa — vl and one transforms from the momentum variables Pa,Pb 
to two new variables one of which is iTa — Vb- This can be done by the 
transformation 

Pia = Wi + -^Ai 

rria 

Pib = Wi--^Ai 
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7oPo{ma + mb)c 

. - - Via- Vib 
= Pia - Pib = n 

7oPoC 

Pia 



Pi, 



Pib 



Pib 

1 1 1 



fi rria rUb 
d%d% = d^Wd^A 



(4) 



Aj is proportional to the relative velocity, iTa — Vb when the velocities are 
non-relativistic. A similar transformation is used in Ref.l and Ref.4 except 
that for them the particles are identical and the transformation is simpler. 

S{piaPja) can be written as 

S{PiaPja) = PiaQja + PjaQia + QiaQja 

Qia = P'ia-Pia (5) 

This result can written as 



1 , r/ 

(lia)- 

Qia = qia/iloPolJ'C) 



S{PiaPja) = [{WiQja + M^j^ia) ] + [( fi^Qja + AjQia + ?ia?ja)] 



(6) 

Eq.3 can be rewritten in terms of W, A as 

<5{PiaPja)> = ^6 / d^xd^Wd^Afa{x,Pa)h{x,Pb)\Va-Vb\ 
GabdO! S{piaPja) dt 
KPiaPja) = [{WiQja + ^j^ia)—] + [{ — f{Aiqja + AjQia + qiaqja)] 

(7) 

One may note that (Tab depends only on A and not on W. In the expression 
for S {piaPja) the second term will be seen to depend only on A and gives rise 
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to the usual intrabeam scattering growth rate, while the first term depends 
on W and will be seen to vanish for identical particles and gives rise to the 
cooling rates for ion electron scattering. 

The transformation from Pa,Pb to W, A allows us to do the integral over 
dQ'. Eq.7 holds in any CS where the particle motion is non-relativistic. For 
each Pa,Pb one can define a center of mass CS, called the CMS, in which 
p^+pt,^ 0. In the CMS 

= Pia - Pib = PialilOpOlJ^c) 

In the CMS, A and pa have the same direction, and pa is scattered by the 
electrons to Pa ' which is along the direction given by the polar angles 6*, 
relative to the direction of Pa or A. 

In Eq.7, only the depend on the scattering angles ^,0 .To do the 
integral over dQ,' in the Rest CS one has to evaluate the integrals 

di 

Cij 

dfl'aab is an invariant and A, (ja are both the same in the CMS and the Rest 
CS as they are both the difference of 2 vectors that are proportional to a 
velocity, di, Cij are tensors in 3-space. If these integrals are evaluated in the 
CMS and the result is written in terms of tensors in 3-space then the result 
will also hold in the Rest CS. 

In the CMS, we introduce a polar coordinate system 9, where 9 is mea- 
sured relative to the direction of or A and we assume that aab{9, (p) is a 
fumction of 9 only, we can then write 

A = (0,0,1)|A| 

p-a = (0,0,l)|A|(7o/3o/ic) 

Pa' — (sin COS 0, sin^^ sin 0, cos^^)| A|(7o/5o/ic) 

Qa — (sin6'cos0, sin6'sin0,cos6' — l)|A|(7o/3o/xc) (9) 

In the CMS, using Eq.9, one finds 

di = -27r [ d9sin9{l-cos9)aab{0,0,l)\^\ 



J dfl'aabQia 

j dn'aabli^iQja + ^j^ia) + QiaQja] 



(8) 
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Cij 



ttJ^ dOsiv^eaah 1 




(10) 

In computing Cij one may note that the Ajg^a + ^jQia term in Eq.8 only con- 
tributes to C33 while the qiaQja term contributes to to all 3 diagonal elements 

of Cij. 

These results for d^, Cij in the CMS can be rewritten in terms of tensors 
in 3-space as 



di — —2tt j d0sin9{l — cosO)(Tab^i 

■■ d9 sm^ 9aab ( | A| ^5^^ - 3Ai ) 



Cij 



(11) 



In this form the results will also hold in the Rest CS. Eq. 7 can now be 
rewritten as 

<S(piaPja)> ^ ^bj d^xd'^Wd'^Afa(x,Pa)fb{x,Pb)\Va-Vb\ 

{[-27r—{WiAj + WjAi) f desine{l - cose)aab]i 
+[7r(-^)'(|A|'(5y - 3AiA,) / desiv?eaab\2) dt 

(12) 

Eq.l2 can be used to compute either intrabeam scattering for identical 
particles or electron cooling. If the a and b particles are identical, then the 
second term indicated by [ ]2 and called the A-term gives the growth rates 
for intrabeam scattering. In this case, the first term, indicated by [ ]i and 
called the W-term, will vanish. This is shown below for gaussian distributions 
and also can be shown to hold for any distribution because of the symmetry 
of the a and b particles. If the b particle is much lighter than the a particle, 
the W-term gives rise to coohng of the a particles and the A-term is smaller 
than the W-term by the factor irib/ma. This is shown below for gaussian 
distributions. Eq. 12 holds for any distibutions, fa{x,Pa),fb{x,Pb)- In the 
next section, we will specialize to gaussian distributions. 
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it is often assumed that aab is given by the Coulonb cross-section in the 
CMS CS for the a and b particles. This is given by 



(^ab = (^j 



rab 



ZaZbe^ 



PabC = \Va-Vb\ 

(13) 

The integrals over 9 in Eq.l2 can then be written as 



/ 



dOsinOil — cos9)-. -r- 

(1 — cosOy 

1 



In 



I 



1 - cosef 
tan{dmva) 



1 + 

In 
rab 



'(^Ibk 



'ab^'^naxab \ 

rab ) 

l^ab^Tiiaxab \ 
. ^-ab J 



1 + {rab/iPlh^ 



maxab 



ab 'maxab 



bmaxab IS the maximun allowed impact parameter in the CMS. 9min is the 
smallest allowed scattering angle in the CMS. 

It will be seen below that to compute the cooling rates for the emittances 
one will also need the cooling rates for < XiaPja >■ When the a and b particles 
are identical, the < XiaPja > are zero , but not zero when the particles are 
different. Using Eq.7, one finds 

<S{XiPja)> ^ Nbj d^xd^Wd^Afa{x,Pa)fb{x,Pb)\Va-Vb\ 

UabdO! S{XiaPja) dt 

r- _ - 1^ 
X^OPja — XiQja' 



rUn 



5{XiPja) 

From Eq.ll one has 

j dVt'aabQja — ~27r j d9sin9{l — cos9)aabAj 
which gives 



(14) 
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<S{XiPja)> = Nb j d^xd^Wd^A fa{x,Pa)fb{x,Pb)\Va-Vb\ 

[— 27r j d9sin9{l — cos9)aab XiAj^^] dt 

(15) 

Eq.l5 shows that < S{xiPja) > gives rise to a coohng term which vanishes 
when the particles are identical , or when ai — for the ion partcle for a 
gaussian distribution. 



4 Cooling rates for < PiaPja > in the Rest CS 
for Gaussian distributions 

In this section, we will find the cooling rates due to the scattering of the ions 
by the electrons in the cooling section when the ion and electron bunches 
have gaussian distributions. In Eq.l2 , we will keep only the 1^-term as the 
A— term , discussed later, is smaller by the factor nib/ ma In this paper, it 
will be assumed that the dispersion is zero in the cooling section. 

For a gaussian distribution, fa{x,Pa) is given for the ion bunch for zero 
dispersion by Ref.[3], 

faix,Pa) = ^exp[-Sa{x,Pa)] 



j d^xd^p exp[-Sa{x,Pa)] 



^xa ~l~ ^ya ~l~ ^sa 
1 



(16) 



'^xai^i X^) Xg^ — Pxa/ POa 
^xa 

li^, Xg^ [x + (^Px-^a "l~ Q^ia-^) ]/ Pxa 



Sya = —eya{y,ya) v'a = Pya/POa 
^ya{y,ya) = + Wyya + Oiyayf]/Pya 
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Ss 

es{s,Ps/pOa) 

^s{S:Ps/POa) 

es{s,Ps/pOa) 
Ps 



—es{s,Ps/poa) 

^ (Ps/POa) 



2(7? 



Ps 



[S^ + PsiPs/POa) 



= [{sf + mPs/poam/Ps 



(17) 



A longitudinal emittance has been introduced so that the longitudinal motion 
and the transverse motions can be treated in a similar manner. (5s in the 
Rest CS is larger than 13s in the Laboratory CS by the factor 7q. s,Ps are 
the paricle longitudinal position and momentum in the Rest CS. 

In Eq.l2 we will now do the integration over d^xd^W using the above 
gaussian ditributions. Because there is no dispersion in the cooling section 
the integral over dxdW^ or dsdWs or dydWy can each be treated in a similar 
way. Eq.l2 can now be written using the Coulomb cross-section as 



S < {PiaPja) > = 



I d^xd^Wd^Aexp[-{Sa + Sb)] \Va - Vb\ 

I "''ah N 2 



In 



rria ' "P; 
-27T{WiAj + WjAi) 



1 + 



'Plbb: 



niaxab 



rab 



dt 



(18) 



We rewrite Sa + Sb as 

Sa + Sb ■ 
Sin 



'^i{Sia + Sib) i — X-iVi S 
2 



1 



Sin 



1 



P 



Pia rUa 



P 



1/2 



Sia + Sib = Aiiixf + A22iW^ + 2AuiXiWi + {A^iXi + AouWi)Ai + ^ooi A: 
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Aqu — 



m ti 



A 



22i 



m ti 



(19) 



The symbols [( )]+ and [( )]_ are defined by 

[( )]+ = ( )a+( )b 
[( )]- = ( 

Wc will now make a transformation to eliminate the 2Ai2iXiWi term in Sia + 
Sib- We rewrite S'ja + Sib ^ 



ib 



[Aux^ + A22W'' + 2Ai2xW + (^lox + Aoiiy) A + ^ooA']i 

^12 



^12- 



[^'(^11-X^) + (^2^V + -3^X)^ 
^22 

+(^iox + Aoiiy) A + AoA'li 



tX 



Prji 



A12 

22 T ^-^^2^ 



+ 



Xi 



dxidwi 



A^2[i 



X. 



7)1 



4 1/2 
^22 - 

[^11^2 

■ 4l/2' 
^22 

Al/2 



22 









1 


[ AV2J 


i 


.1/2 

L^22 J 



Sia + 'S'ift 



[77^ + + (Aiox + AoiW^) A + AooA% 
[rj^ +pI + {Biorj + BoiPr,)A + AooA% 
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Bioi — 



A.22 . Ai2 







-Boli — 

i 


L ^22 J 



(20) 



In the expression for Sia + Sn, given at the end of Eq.l9, the hnear terms 
in T], can be ehminated by the transformation 



Prii 



, -^01 A 



Sia + Si, = +pV + (^o-Sio/4-So\/4)A2], 



Eq.l7 can now be rewritten as 

m 1 
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A)],A, + [W^P.(P.-^A)],A, 



(21) 



< S{piaPja) > 



J d^fid%d^A exp[-(Sa + Sb)]\Va-Vb\ 

dt 



maxab \ 



^1/2 
P 



'"'a Pab 

-27r[[W,{v - ^A)],A, + [W,^{p-, - ^A)],A, 
+[Wr,{fj - ^A)]jAi + [W^^ipr, - ^A)],A,] 

^1/2^1/2^1/2 



(22) 



Using Eq.20 for Sia + Sn, and for Wij , one can do the integral over d^rjd^Pri 
and get 
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TaTbA'i' ma 



"J 



/ 



In 



1 + 



maxab \ 



Tab 



Boi 



Pv 



dt 



= 7o/3o(A^ + A^ + A^)V^ 



B 



00 



01 \2 



[AiiA 



22 



A' 



11/2 41/2 11/2 



X 



Bioi 



A 



lOli 



A 



1/2' 
22 



ryi 



Al/2 



A 



12 



Al/2 

-Soil - 



Pr)i 



4 1/2 
^22 J 



[A^iW, 



A 



4 1/2 
^22 



io]ji7i 



12 



01 



^V2 



-Boli — 



iOl 



4 1/2 
^22 



A 



22i 



Al2i = 




^lOi — 


2 At Oi" 






+ 


. m Ci . 



m e,; 



r A^\2 A 



-I + 



(23) 



Eq.23 is our final result for the cooling rates for < PiaPja > in the Rest 
CS, for two overlapping gaussian bunches , with no dispersion in the cooling 
section. For this case one gets zero results when i ^ j. The remaining 3- 
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dimensional integral over d^A is an integral over the relative velocities of the 
ions and electrons. 

It will be seen below that to compute the cooling rates for the cmittances 
one will also need the cooling rates for < XiaPja >■ For gaussian distribu- 
tioins, using the coulomb cross section and Eq.l5, Eq.l8 is replaced by 

S<{XiPja)> = :z— ^ / d^xd^Wd^Aexp[-{Sa + Sb)]\Va-Vb\ 

i- ni- h J 



-X 



Tab 



H 2 



ah . 



In 



1 + 



'maxab 



Tab 



dt 



X 



-2nXiAj 



(24) 



After going from the x, W coordinates to r],Pn and integrating over 77, 
Eq.23 is replaced by 



5 < {XiPja) > 



X 



Nb 1 3 2 1^ ^ 

rTT r„^C Xi 



J ^3^ exp[-{X,Al + XyAl + XsAj)] ^ ^ 



In 



1 + 



/32 t \ ■ 

H ab "unaxah \ 



Tab J 



dt 



X. 



BlOi 



A. 



1/2' 
22 



iV2 
[AioXr^ + AoiWr, 



A. 



1/2 
22 



iV2 



-A 



A 



A 



12 



iV2 



Prfi 



4V2 
^22 



OH 



12 



01 



AV2 



[^01 W^, 



Pi) J 



oil 



^01-172 

^22 



(25) 



5 Emittance growth rates 

One can compute growth rates for the average emittances, < > in the 
Laboratory Coordinate System, from the growth rates for < PiaPja > in the 
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Rest Coordinate System. In the following , dt is the time interval in the 
Laboratory System and dt is the time interval in the Rest System, dt — ^dt. 
The final results are, for zero dispersion, 



d _ Pia d -2 , '^(^ia d 



(26) 



To derive the above results, the simplest case to treat is that of the vertical 
emittance. The verical emmitance is given by 

^ya{y, y'a) = [y'^ + {Ma + ayayf]/l3ya y'a = Pya 

Seya = (3yaS{plJ + 6{2ayay{pya) 
dt ^''"^ 7 di^^'^ ^ ^ 7 dt y^'" ^ 

(27) 

In Eq.(27), = pya, Seya is the change in Cya in a scattering event. Similar 
results will hold for exa and e^a for zero dispersion. 



The A term in electron cooling 

In the previous section it was assumed that in Eq.l2 one could drop the 

second term or A term compared to the first term or W term. This is true 
when nib << fT^a and pa — pb in the Rest CS. Using Eq.4, one can write 



A,; 



. ma + rrib 



+ Pb 



nib 



ma + nib 



= \pa- Pb]i 

(28) 

Thus W and A are both of the same order as pa ■ If the motion is non- 
relativistic in the Rest CS, ~ A ~ pa- From this it follows that the A 
term in Eq.l2 is smaller than the W term by the factor mb/nia- 
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It has also been assumed that the motion in the Rest CS is non-relativistic. 
In the Laboratory CS, the rms spread in the relative momentum is given by 



'pi 



1/2 



x,y,s 



(29) 



For gold ions in RHIC at 7 = 100 



= = 5e — 8, (3x = I3y = 50 and = CTpy = 2.24e — 5 

eg = 1.8e — 4, Ps = 300m and aps = .55e — 3 

In the Rest CS cXpa;, apy are unchanged at 2.24e-5 And aps is reduced by the 
factor 7 to .55e-5 The spread in each of the momenta in the Rest CS is of 
the order of le — SrUaC since 7 = 100 and the ion velocities are of the order 
of le-3c. Similar numbers hold for the electrons in the electron bunch. 
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